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Abstract 

We calculate the a s corrections to the transverse polarization fraction of ?/>"s 
produced at the Tevatron. If the '^'-anomaly' is explained by gluon fragmen- 
tation into a colour octet Fock component of the ip', the ip n s should be 100% 
transversely polarized at leading order in a s , up to spin symmetry breaking long- 
distance corrections. We find that the short-distance correction to the transverse 
polarization fraction is a few percent, so that a polarization measurement would 
provide a reliable test of the colour octet mechanism. 
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1. Introduction 

The topic of quarkonium production has received renewed attention as of late due to a 
plethora of data and new theoretical developments. While it has been known for quite some 
time that quarkonia calculations are tractable as a consequence of the fact that the non- 
perturbative dynamics may be swept into non-relativistic wave functions, an organizational 
principle was still lacking. This shortcoming manifested itself in the fact that the theory of P 
wave decays was plagued by infrared divergences that could not be factorized. This problem 
was resolved KJ once the calculation was organized using non-relativistic QCD (NRQCD) 
@. Within the framework of this effective field theory, the scales involved in quarkonia 
calculations are cleanly separated and factorization becomes explicit. 

NRQCD has been widely used in calculations of quarkonium production rates. For trans- 
verse momentum of the quarkonium larger than the quarkonium mass, the dominant 
production mechanism will be fragmentation. This is essentially a resonance effect resulting 
from the fact that the off-shellness of the fragmenting parton is small (~ 4m;?) compared 
to the transverse momentum squared. The rate for this process is enhanced by a factor of 
p 2 -, / (4m 2 ) , which can overcome the suppression due to any additional powers of a s present 
in the fragmentation amplitudes. The fragmentation functions for quarkonia are partially 
calculable in perturbation theory ||. 

If one does not include the fragmentation contribution, the prediction for the prompt 
production cross sections for both J/ip and ip 1 at the Tevatron are off by several orders of 
magnitude. The inclusion of fragmentation into color singlet quark-antiquark pairs brings 
the J/ip production rate into agreement with experiment within theoretical errors j|, [5|, [|. 
but still leaves the theory an order of magnitude short for ip 1 production ||. A possible 
explanation |7| for this discrepancy is that the dominant contribution to ip' production 
results from fragmentation into two charmed quarks in a relative color octet state.Q This 
is easily seen once we consider the process within NRQCD. In this formalism the full Fock 
state decomposition of the if}', 

W) = |ccflsf >]) + 0(v) \cc[ 3 Pi%) + 0(v 2 ) + 0(v 2 ) IcefS^jgg) + ..., (1) 

is taken into account by the set of all possible production matrix elements [K|. The nota- 

1 Recent measurement of the J/ip production fraction from \c decays indicates that colour octet production 
is also important to account for the fraction of J/ip's not coming from B and \c decays [|[ ^. 
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tion here is spectroscopic 2S+l Lj, with an additional superscript in parenthesis to denote the 
colour state. Each Fock component is weighted by some power of the heavy quark veloc- 
ity. This weight is determined by the velocity scaling rules which will be discussed below. 
Although the transition of a cc state in a relative colour octet state into a ip' is suppressed 
by v 4 ~ (0.3) 2 relative to the transition from a colour singlet pair, this suppression is off-set 
by an (7r/a s (2m c )) 2 -enhancement for the production rate of cc pairs in a 3 S[ 8) state. Ad- 
ditional numerical factors ensure that the value of the matrix element (associated with the 
hadronization of the colour octet pair) required for agreement with experiment is indeed 
consistent with the velocity scaling rules J7J. 

As noted by Cho and Wise > the colour octet production mechanism could be checked 
by measuring the polarization of the ip' , as they should be strongly transversely polarized at 
leading order in a s . The reason for this asymmetry is that the fragmenting gluons are nearly 
on-shell, so that the cc pair inherits the gluon's transverse polarization up to corrections of 
order (4m 2 )/g 2 , where go is the gluon energy in the lab frame. Due to spin symmetry of 
the leading order NRQCD Lagrangian, the polarization of the quark pair stays intact in the 
subsequent (nonperturbative) evolution into a ip' state up to spin symmetry breaking correc- 
tions of order v 4 /3 ~ (3 — 4)%. In this paper we address the question whether this signature 
persists after inclusion of radiative corrections to the short- distance production process and 
calculate the order a s corrections to the prediction of 100% polarization (neglecting spin 
symmetry breaking). We find that these corrections are small, and thus the colour octet 
proposal can be easily tested since the polarization will not be destroyed by radiative cor- 
rections. Other deviations from pure transverse polarization from (spin-symmetry breaking) 
long-distance effects or higher twist effects will most likely be smaller than the short-distance 
correction, except, possibly, at small pr, where higher twist effects can become large. 



2. Fragmentation functions 

The fragmentation of a gluon (or any parton) into a quarkonium is a two-scale process. 
Over distances of order l/m c the gluon fragments into a cc pair with longitudinal momentum 
fraction z and small relative velocity v in the rest frame of the pair. The subsequent arrange- 
ment of the quark pair into a quarkonium bound state takes place over longer distances, of 
order l/{m c v). To the extent that « 2 C 1, the effects on both scales can be separated in 
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NRQCD and the longitudinal momentum fraction of the quarkonium identified with that of 
the quark pair, so that the ^-dependence is perturbatively calculable. The fragmentation 
function, defined for instance as in fl2| , can then be written as 

D g ^ Hx (z,v) = (O^W (2) 

n 

where H\ denotes the quarkonium in a specific polarization state A. The d n -ij are short- 
distance coefficients, independent of the quarkonium state and computable as expansions 
in a s . All bound state information is summarized in the vacuum matrix elements of the 
operators 0^.h. These parameters must be determined phenomeno logically or from lattice 



calculations. Following the notation of ||10|| , the generic form for the production operators is 

Onk = X^n-,i i/> a ( H ] a { H ] ift Kf n . d X, (3) 

where if), \ are nonrelativistic quark and antiquark fields and destroys a quarkonium state 
H with polarization A in the out-state. The matrix elements are supposed to be evaluated 
in the quarkonium rest frame. The kernels n' n .- contain colour matrices, spin matrices 
and derivatives. The colour and spin indices are suppressed. The operators O u \a depend on 
two sets of three- vector indices i and j. If we consider unpolarized production and include 
summation over A in the definition (||), the matrix elements (O^.^) can not depend on any 
three-vector. A tensor decomposition then reduces the sum in (Q) to a sum of a few scalar 
operators as listed in flOfl . In the case of polarization this decomposition also includes the 
polarization tensor of the quarkonium. After decoupling of indices from the short-distance 
amplitude, the kernels K n contain projections on the spin and angular momentum orientation 
of the quark pair. The dependence on the factorization scale for NRQCD, m c v < fx' < m c , 
drops out in the sum (Q), while the dependence on /i, which separates the gluon production 
from the fragmentation process is governed by the usual Altarelli-Parisi evolution. 



3. Colour octet production at leading order 

Given the Fock state decomposition ([]]), the leading contribution in v 2 comes from gluon 

and charm fragmentation into the colour singlet 3 5i state. Charm fragmentation is shownQ 

2 This diagram as well as all others are understood to be calculated in axial gauge. Otherwise, diagrams 
with gluons attached to eikonal lines would also have to be considered, see |12]. 
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in Fig. p]b. Everything prior to the shaded oval is considered as short- distance dominated, 
whatever is to the right is contained in the matrix element of an operator . In Fig. [I]b, 
this is simply the non-relativistic wave-function (squared) at the origin. Using the velocity 
scaling rules derived in colour singlet charm and gluon fragmentation yield contributions 
of order a 2 v 3 and a 3 t> 3 respectively. The net contribution from singlet fragmentation falls 
short of the experimental data by a factor of about 30 [||. However, since v 2 ~ 0.3 is 
not small, the contribution from operators, whose matrix elements are suppressed in v 2 , 
may be important. This is especially so for 0$ ( 3 Si), because it is the only operator with 
a short- distance coefficient d(z) at order a s . A quark pair in a colour octet 3 S\ state has 
non-zero overlap with a Fock state component of the ip' that contains two additional gluons, 
represented by the emission from the black boxes in Fig. |l|a. The long-distance dynamics is 
accurately described by the NRQCD Lagrangian and the velocity scaling rules apply. The 
quark pair in the 3 S{ 8) state can evolve into the ip' by a double electric dipole transition, each 
of which contributes a factor off to the amplitude as indicated in (|1|). Thus, the contribution 
from colour octet gluon fragmentation starts at order a s v 7 . The short-distance coefficient 



was computed in |TJj] to leading order in a s and the fragmentation function is given by 

= " *) (Of(%)) (4) 

with (OffSi)) defined as in [[IIJ. Braaten and Fleming found |7| that the prompt if)' 
production cross section at the Tevatron and its pr-dependence can be fitted by (Of ( 3 Si)) = 
0.0042 GeV 3 , a value consistent with a v 4 -suppression relative to (Of ( 3 Si)) = 0.11 GeV 3 . 

Since the 1/p^-dependence of the production cross section is rather generic for any frag- 
mentation mechanism, additional information is necessary to establish the colour octet mech- 
anism. One important signature of the colour octet mechanism is the ip' polarization. The 
polarization of the quarkonium is calculable to a certain degree, because the non-relativistic 
dynamics that governs the long-distance evolution conserves the heavy quark spin up to cor- 
rections of order v 2 that enter through the magnetic coupling a-B in the NRQCD Lagrangian. 
To implement the implications of heavy quark spin symmetry on the tensor decomposition 
of matrix elements, it is useful to represent the quarkonium state as a product of a spin and 
orbital angular momentum part, similar to the decomposition in heavy and light degrees of 
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(b) 






Figure 1: (a) Leading order gluon fragmentation into ccpS^ ]. (b) Charm fragmentation, 
(c), (d) Real and virtual fragmentation diagrams at order a 2 s . The black boxes denote 
nonperturbative evolution of the cc-pair into a if)' through soft gluon emission. 



(5) 



freedom for heavy-light mesons @. For the 3 S[ 8) -operators one obtains 

i(Or( 3 5 1 ))e a *(A)e fe (A) 

{Of ("Si)) Ex e«*(\')e b (\') |(1(A - A'); 1A'|JA>| 2 
where e a (A) is a spin-1 polarization vector in the quarkonium rest frame. The first line holds 
for H = ip' and the second for H = Xcj- When inserted in (Q), these equations determine 
the required projections of the short-distance amplitude on the spin (and, in general, orbital 
angular momentum) of the heavy quark pair with appropriate Clebsch-Gordon coefficients. 
They are equivalent to and reproduce the spin symmetry constraints investigated in [|TT| . 
For a ip' with polarization A, one must project on a quark pair with equal polarization, as 
expected, since the two El transitions do not change the spin. Since the quark pair in Fig. [TJa 
is transversely polarized, the projection with e a (0) vanishes, and, at leading order in a s , we 
get 



D g _M 0, n) = Dg-*^. (z, y) = Dg^> {z, fx) 



(6) 
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for the fragmentation functions into longitudinally and transversely polarized ^"s. Because 
spin-symmetry is not exact, long-distance corrections of order v 4 arise, if the quark pair in 
the octet state undergoes a double magnetic dipole transition. 



4. Short-distance corrections to transverse polarization 

Corrections to pure transverse polarization also arise, when the polarization of the quark 
pair prior to nonperturbative evolution is modified by the emission of hard gluons. The 
corresponding real and (some of the) virtual corrections to gluon fragmentation are shown 
in Fig. 0c and d. Let us define 

« sl -rTZ sl -« (7) 

as the ratio of the (prompt) production cross section for transversely and unpolarized ip n s. 
Since L vanishes at leading order, 5£ is obtained from the fragmentation function into longi- 
tudinally polarized tp n s at order a^, while the total production rate L + T in the denominator 
can be evaluated in leading order approximation. 

For the same reason as in leading order, the virtual corrections do not yield longitudi- 
nally polarized ip n s. This results in considerable simplification, since along with the virtual 
corrections, radiative corrections to the matrix elements in NRQCD would also have to be 
calculated. The full matching calculation for the next-to-leading short- distance coefficient of 
(Of ( 3 5*i)) has been performed by Ma |16| for unpolarized quarkonia. When combined with 
the result for the fragmentation function into ip' L in eq. (|H]) below, the production rates for 
longitudinally and transversely polarized ^"s through a quark pair in a state could be 
obtained separately at next-to-leading order in a s . Here we consider only the ratio £, so that 
the total rate can be evaluated in leading order. 

After hard gluon emission the quark pair that enters the shaded ovals in Fig. [TJc can be in 
various spin and angular momentum states. At order a^v 7 , the three distinct possibilities are 
3 Sf 8) , 3 Pf and ^ 8) . The power v 7 arises as combination of factors of v associated with the 
probability for binding a quark pair in a certain angular momentum state (v 3 for S'-waves 
and v 5 for P-waves) and additional factors associated with the amplitude of the relevant 
Fock state in (§. For 3 S{ 8) and x sf } a factor v 4 is obtained in the amplitude squared for the 
double El and single magnetic dipole transition, respectively. In case of 3 Pj 8 ^ a single dipole 
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transition provides v 2 . 

Since there is no interference between S'-wave and P-wave amplitudes and spin-0 and 
spin-1 amplitudes, the contributions from the three different intermediate quark pair states 
can be determined separately. Let us turn first to the 3 S[ 8 ^ intermediate state. This is the only 
one to which diagram (ciii) contributes. The projections onto the short-distance amplitude 
are determined by (|5|). In the actual calculation it is useful to boost into a frame where the 
'+- component' of the quarkonium momentum is large and to use the covariant projection 
operators of \ 17\ . We then follow the method of |3|, [14| for computing fragmentation functions. 
The result for the integral of the fragmentation function can be expressed as 

i ^ ds 1 

J dz d(z, /i, //) = J — J dzA(s,z), (8) 

with A(s, z) a certain projection of the amplitude squared. After interchange of integrations, 
one can read off the fragmentation function. In general, the above expression is divergent. 
Divergences come from the upper limit of integration over invariant mass s of the fragmenting 
gluon (cut off by the fragmentation scale /z above) due to the collinear approximation implicit 
in the definition of the fragmentation function. Furthermore, the integral over z may diverge 
due to soft gluon emission. 

For the unpolarized ip' fragmentation function the logarithmically divergent part of the 
integration over s gives 

D%{z^^>) = ^ (Of(%)) • g^(*)ln£^ (9) 

where Pg^ gg (z) is the Altarelli-Parisi gluon-gluon splitting function. The logarithmic diver- 
gence occurs as part of the production amplitude, when the fragmenting gluon splits into 
two collinear gluons of low invariant mass. Since both gluons are parallel and transverse, so 
is the final quark pair. Thus, the s- integral is finite for longitudinally polarized r/>"s. 

For a quark pair in a 3 S[ 8 ^ state the integral over z is finite. A divergence can poten- 
tially arise for P-waves, when the amplitude is expanded in the relative velocity. However, 
soft gluon emission with energy less than //, the NRQCD factorization scale, should be 
considered as part of a NRQCD matrix element. For longitudinal polarization there is no 
contribution at order a s , into which an infrared logarithm could be absorbed and all longi- 
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tudinal fragmentation functions must be finite as z — > 1. The integral above can then be 
evaluated taking fi to infinity and // to zero in all cases considered in this paper. 

The result for the fragmentation function into longitudinally polarized t/>"s through a 
colour octet 3 S\ quark pair is 

fl * (A) ^^tfM, (10) 

y Vl 24 z m A c 

where iV c = 3 is the number of colours. 

For a ip' produced through an intermediate quark pair with orbital angular momentum 
one, we need the matrix elements 

( X ^ a T A f-j Di) i> 4 A)t «t A) ^ b T A (-- D 3 ) x) = (Ot\ 3 P )) Su e a *(X)e\X). (11) 



v 2 "7 r r " V 2 / 

The right hand side incorporates the constraints from spin symmetry and is accurate up to 
f 2 -corrections. Thus, we have to project the short- distance amplitude on a longitudinally 
polarized quark pair and sum over all orbital angular momentum states. Note that the total 
angular momentum J of the quark pair is not specified. The resulting fragmentation function 
is 



z 2 + 8z - 12 



' w 16 iV c 

2z 3 + z 2 - 28z + 24 



z 2 



1 - z) ln(l - z) 



<°ftf°». (12) 



mi 



2z 

Because of the decomposition of the matrix element (|TTD, which would differ for Xcj states, 
this fragmentation function can not be obtained as a linear combination of the short-distance 



part of the fragmentation functions into polarized Xcj states [18]. Note also that it is finite 
ctS Z ^ X . 3jS expected from the considerations above. 

Finally we have to account for the quark pair in a spin-0 state. While in the previous two 
cases, the nonperturbative evolution could be considered as spin-conserving at order v 7 , a 
non-zero value for (Of (%)) arises through a magnetic dipole transition at order v 7 . Since 
the initial quark pair evolves into a ip' with a probability independent of the polarization of 
ip' , the fragmentation function for a longitudinally polarized ip' is given by dividing the total 
fragmentation probability by three. The result is 

D «?> u ~, 2mA = «'w< : 4 [ 3 , _ 2 - + 2d - ,> mi - 2 ,i MM. (13) 

J m s c 



a-*i> L ^ ' c! 864 N r 
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The functional dependence is the same as for the colour singlet fragmentation function for 
Vc in [§. 

The fragmentation function D g ^^ L (z) for longitudinally polarized ^"s is given by the 
sum of (|10|), (|l~2]) and (|i~3| ) at leading non- vanishing order in a s . There is also a contribution 
from charm fragmentation, which is formally of the same order a 2 s v 7 , where the quark pair in 
Fig. [l|b is in an octet state. The corresponding fragmentation functions are identical up to a 
colour factor to the colour singlet charm fragmentation functions, which contribute at order 
a^v 3 . Since colour singlet fragmentation (including charm fragmentation) can account only 
for 1/30 of all ^"s, colour octet charm fragmentation contributes at most a tiny fraction to 
the total production rate. It can be neglected without affecting our conclusions. 

5. Results 

The production cross section for longitudinally polarized ip n s at the Tevatron is obtained 
by convoluting the fragmentation function with the gluon production cross section in pp 
collisions with center-of mass-energy ^ = 1.8 TeV, 

d 1 d 

JL(pp^^ L + X )= J d z^-(pp^g(P T /z)+X)D g ^ L (z^). (14) 

The longitudinally polarized production fraction <5£ is given by dividing this expression by 
a similar expression with D g ^> L (z) replaced by D g ^(z) given in (||). We have used the 
CTEQ3L parton distributions to evaluate the gluon production cross section. The dominant 
production channel is gluon-gluon fusion, while the qq channel is negligible. The gluon 
fragmentation functions are evolved from the initial scale 2m c = 3 GeV to fi of order pt by 
the usual Altarelli-Parisi equations. To very good approximation, mixing can be neglected 
in solving the AP equations. As input we have chosen a s (2m c ) = 0.27. In addition we have 
implemented the pseudo-rapidity cut \r)\ > 0.6 as in the latest CDF analysis 

The final result for <5£ can be represented as the sum of the three contributions to longi- 
tudinal fragmentation discussed earlier: 

*£»(!*, = r l{PT ,») +r2(pT ,,)-lgjgl + r 3 (p„.)||g|. (15) 
We have added a subscript 'SD' to indicate that spin symmetry breaking long distance 
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001 
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0.001 


0.047 
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0.018 


0.025 


0.001 


0.043 


18 


0.017 


0.024 


0.001 


0.041 


21 


0.015 


0.024 


0.001 


0.039 


24 


0.015 


0.023 


0.001 


0.038 



Table 1: The coefficients as a function of px with fi = Pt and |r/| > 0.6. The 
last column gives the deviation from pure transverse polarization, 5^sd = ?"i + ?"2j for 
(Of ( 3 Po))/(C , s'( 3 5i)) = 1 and neglecting the contribution from (Of (%))■ 



corrections of order a s v 11 are not included here. The coefficients are listed in Table [I] for 
fi = Pr- 

As can be seen from the Table, the contribution from r 3 is small and can be dropped 
for any reasonably expected ratio of the matrix elements that multiplies r$. This leaves us 
with one ratio of matrix elements which is not very well constrained at present. The velocity 
scaling rules tell us that this ratio should be of order unity. A slightly more definite, but 
crude estimate can be obtained by the following argument borrowed from [0J : The operators 
Of ( 3 P ) and Of( 3 S 1) mix under renormalization [[L(J. The solution to the renormalization 
group equation in the one-loop approximation is given by 



(Of (%))(»') 



(ofCs 1 )M) + 6 ^- 4) - aM) 



In 



(Of( 3 Po)), 



N c (3 m 2 c a si/j,') 



(16) 



where we used (Of ( 3 Pj)} = (2 J + 1) (O| ( 3 P )), as follows from (JTJ), and (3 = UN c /6 - 



Nf/3. The second term is formally enhanced by a logarithm if /j,' and pi are very different. 
Neglecting the first term on the right hand side and choosing the initial scale of order of the 
typical non-relativistic momenta fi' ~ m c v and // ~ m c , we get, with a s (1.5GeV) = 0.35, 
the estimate 

^.la. (i7) 

Although the logarithm is in fact not large, a comparable estimate for \ c states as in ]I[ 
is only off by a factor of two from the phenomenologically known ratio in that case. If we 
assume that this difference arises from an inaccurate choice of the initial scale and adjust the 
logarithm to the known ratio for matrix elements in the Xc case, the estimate above changes 
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Figure 2: The short-distance correction 5£,sd as function of p?- The solid line shows the 
prediction for // = p T and (Of ( 3 P )) / (Of ( 3 Si)) = 1. The two dash-dotted lines are for 
\i = Pt/2 (upper curve) and \i = 2p T (lower curve) at fixed (Of ( 3 P )) / (Of ( 3 5i)) = 1. The 
two dashed lines correspond to variation of (Of ( 3 Po)) / (Of ( 3 Si)) between 0.5 and 2 at fixed 
H = Pt- 

from 1.25 to 0.65. In view of this we consider a ratio of one as our best guess. With this 
choice the prediction for S^sd is shown in the last column of the Table and as solid line in 
Fig. |2|. The Figure also shows the prediction if we modify the above ratio of matrix elements 
to 2 or to 1/2. This uncertainty constitutes the dominant source of theoretical error until 
(Of ( 3 Po)} is determined from another process. For comparison, we have plotted the varia- 
tion of 5£, when the renormalization scale is set to \x = Pt/2 or /x = 2pT- This uncertainty 
is about ±20%, except at small px, where the choice pr/2 is not very well motivated. It 
should also be mentioned that by nature of the fragmentation approximation, the prediction 
is more reliable at large Pt- If one integrates the production cross section over pt starting 
from some pT.min, the corresponding integrated 8£ is roughly equal to <5£ at PT.min, because 
the absolute production cross section decreases rapidly with pt- 

6. Discussion 

From the above we conclude that short-distance corrections to the fragmentation func- 
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tions decrease the transverse polarization of 100% at leading order in a s by approximately 
5% and hardly more than 10%, even if the 3 Po-matrix element is large. We have ignored 
the colour singlet production mechanisms, since they account for only 3% of all i[)"s. If we 
assume that the singlet mechanisms produce randomly polarized ip h s, they would contribute 
only an additional percent to the longitudinal polarization fraction. A potential correction of 
order Am 2 c /p^ comes from the fragmentation approximation. However, the octet production 
processes at leading order in a s yield basically pure transverse polarization even without 
the fragmentation approximation to the matrix elements due to a particular structure of 
the gluon-gluon fusion amplitude ||. Thus the higher twist correction is suppressed by 
(a s /it) 4m 2 /p 2 ^, which is small for sufficiently large px ~ 10 GeV. 

Deviations from pure transverse polarization arise from spin symmetry breaking already 
at leading order in a s , when the quark pair in Fig. |Tja undergoes spin-changing nonperturba- 
tive evolution. Because of parity and charge conjugation symmetry these corrections scale 
with v A . Since random polarization would give <5£ = 1/3, a natural estimate of the long- 
distance correction is 5£ld = 0(v 4 /3) ~ (3 — 4)%. Despite their smallness, the short distance 
corrections constitute the dominant source of depolarization. We estimate, conservatively, 
that if the colour octet mechanism does account for the observed prompt if)' production cross 
section, then the if) h s should be (85 — 95)% transversely polarized. 

This prediction can be tested rather simply by measuring the lepton angular distribution 
in the if)' — ► l + l~ decay channel, 

dT 1 -35£ 2 

oc 1 H — cos 6, 18 



d cos 6 1 + 5£ 

where 9 denotes the angle between the lepton three momentum in the if>' rest frame and the 
if)' three momentum in the lab frame. As statistics improves, a measurement of the angular 
distribution will become feasible. Time will tell whether such a measurement, together with 
other signatures of the colour octet mechanism, such as those in e + e~ annihilation [20], pro- 



vides a consistent picture of quarkonium production that includes the relativistic structure 
of the bound state. 

Acknowledgements. We benefitted from conversations with Wolfgang Kilian and Mark 
Wise. M. B. is supported by the Alexander von Humboldt foundation. 



12 



References 

[1] G.T. Bodwin, E. Braaten and G.P. Lepage, Phys. Rev. D46 (1992) R1914, ibid. D46 
(1992) R3703 



[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8 
[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 
[17 
[18 



W.E. Caswell and G.P. Lepage, Phys. Lett. B167 (1986) 437 
E. Braaten and T.C. Yuan, Phys. Rev. Lett. 71 (1993) 1673 
M. Cacciari and M. Greco, Phys. Rev. Lett. 73 (1994) 1586 

E. Braaten, M. Doncheski, S. Fleming and M. Mangano, Phys. Lett. B333 (1994) 548 

D. P Roy and K. Sridhar, Phys. Lett. B339 (1994) 141 

E. Braaten and S. Fleming, Phys. Rev. Lett. 74 (1995) 3327 

P. Cho and A. Leibovich, CalTech preprint CALT-68-1988 |hep-ph/ 9505239" | 



M. Cacciari, M. Greco, M.L. Mangano and A. Petrelli, Phys. Lett. B356 (1995) 553 
G.T. Bodwin, E. Braaten and G.P. Lepage, Phys. Rev. D51 (1995) 1125 
P. Cho and M.B. Wise, Phys. Lett. B346 (1995) 129 
J.C. Collins and D.E. Soper, Nucl. Phys. B194 (1982) 445 

G. P. Lepage, L. Magnea, C. Nakleh, U. Magnea and K. Hornbostel, Phys. Rev. D46 
(1992) 4052 

E. Braaten and T.C. Yuan, Phys. Rev. D50 (1994) 3176 

A. F. Falk, Nucl. Phys. B378 (1992) 79 
J. P. Ma, Nucl. Phys. B447 (1995) 405 

B. Guberina, J.H. Kiihn, R.D. Peccei and R. Ruckl, Nucl. Phys. B174 (1980) 317 
P. Cho, M.B. Wise and S. Trivedi, Phys. Rev. D51 (1995) 2039 



13 



[19] The CDF collaboration, Fermilab preprint Fermilab-Conf-94/136-E, contribution to the 
27th International Conference on High Energy Physics, Glasgow, July 20-27, 1994 



[20] E. Braaten and Y.-Q. Chen, Northwestern University preprint NUHEP-TH-95-9 ||hcp- 



ph/9508373fl ; K. Cheung, W.-Y. Keung and T.C. Yuan, Femilab preprint FERMILAB- 
PUB-95/300-T ||hep-ph/ 9509308 1 



14 



